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We construct a solvable spin chain model of many-body localization (MBL) with a tunable mo¬ 
bility edge. This simple model not only demonstrates analytically the existence of mobility edges 
in interacting one-dimensional (ID) disordered systems, but also allows us to study their physics. 
By establishing a connection between MBL and a quantum central limit theorem (QCLT), we show 
that many-body localization-delocalization transitions can be visualized as tuning a mobility edge 
in the energy spectrum. Since the effective disorder strength for individual eigenstates depends on 
energy density, we identify “energy-resolved disorder strength” as a physical mechanism for the ap¬ 
pearance of mobility edges, and support the universality of this mechanism by arguing its presence 
in a large class of models including the random-field Heisenberg chain. We also construct models 
with multiple mobility edges. All our constructions can be made translationally invariant. 
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In one- and two-dimensional disordered free(-fermion) 
systems, Anderson localization says that all single¬ 
particle states are localized, and thus rules out the ex¬ 
istence of a mobility edge in the energy spectrum sepa¬ 
rating localized and delocalized states [2]. As an active 
area of research, MBL studies the effects of interactions 
added to an Anderson insulator [4, 10, 14, 19-21, 25]. An 
important problem is whether mobility edges can exist in 
ID disordered systems in the presence of interactions. If 
they do exist, what is the physics of mobility edges? 

Recently, numerical progress has been made in favor 
of the existence of mobility edges [16-18], although there 
are still different opinions in the community [7]. To date, 
exact diagonalization (ED) is the only numerical method 
for static properties of quantum many-body systems at fi¬ 
nite energy density, but it is limited to small system sizes 
for the simple reason that the dimension of the Hilbert 
space grows exponentially with the system size. In partic¬ 
ular, the state-of-the-art ED is able to work with Hilbert 
spaces of dimension 705432 or < 20 spin-l/2’s [17], and 
finite-size effects are not always negligible. (The time- 
evolving block decimation algorithm [31] can efficiently 
simulate the dynamics of ID MBL systems [3, 29, 35] due 
to the slow growth of entanglement [1, 13, 27, 32, 33].) 

We construct a “solvable” spin chain model of MBL 
with a tunable mobility edge. Our approach is fully an¬ 
alytical and (almost) rigorous so that there should be no 
confusion about the existence of mobility edges in inter¬ 
acting ID disordered systems. Furthermore, this simple 
model allows us to study the physics of mobility edges. 
By establishing a connection between MBL and QCLT, 
we show that many-body localization-delocalization tran¬ 
sitions can be visualized as tuning a mobility edge in the 
energy spectrum. Although the disorder strength in the 
Hamiltonian is fixed, the effective disorder strength (to 
be defined precisely) for individual eigenstates depends 
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on energy density. This motivates us to identify “energy- 
resolved disorder strength” as a physical mechanism for 
the appearance of mobility edges, and we support the uni¬ 
versality of this mechanism by arguing its presence in a 
large class of models including the random-field Heisen¬ 
berg chain. We also construct models with two mobil¬ 
ity edges such that one region of localized (delocalized) 
states sandwiches two regions of delocalized (localized) 
states. All our constructions can be made translationally 
invariant. 

Physical picture. - Before presenting mathematical de¬ 
tails, we discuss intuitions and essential ideas at a non¬ 
rigorous level. To be specific, we consider the random- 
field Heisenberg chain 

H({hi} ie z ) = X 3 ' 3+i - Mf, (!) 

ie z 

where Si = (Sf,Sf, Sf) is the spin-S' operator at the site 
i, and hfs are independent and identically distributed 
(i.i.d.) uniform random variables on the interval [—h, h\. 
We have two different arguments for the existence of mo¬ 
bility edges in this model. 

Argument 1. —Since the ground-state energy and 
bandwidth of (1) depend on {hi}, we take the union of all 
eigenvalues and eigenstates for all disorder realizations: 

U {{E,m : H{{hi}M) = Em, (2) 

{ hi}e[-h,h]® z 

where H({hi}) is called the parent Hamiltonian of \if>). 
We argue that the states in (2) with larger (in absolute 
value) energy densities are more likely to be localized. 

Intuitively, different disorder realizations have different 
disorder strength, and we roughly quantify the disorder 
strength in each individual disorder realization {h t \ by 
JA \hi\/n, where n is the system size. Admittedly, this 
is not a faithful measure of disorder strength for some 
disorder realizations (e.g., hi = h for Vi € Z), but such 
cases are rare. We observe that the parent Hamiltonians 
for the states in (2) with large (in absolute value) energy 
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densities must have large disorder strength. For example, 
we fix 5 = 1/2 and h = 2.5. (It is shown numerically that 
(1) with 5 = 1/2 has mobility edges for 1.5 < h < 3.5; see 
Figure 1 in [17].) Then, a state with E/n < —1 implies 
Yli \hi\/n > 2 + 2 E 0 « 1.11, where E 0 = 1/4 —ln2 is the 
ground-state energy density of the (homogeneous) spin- 
1/2 anti-ferromagnetic Heisenberg chain. 

Argument 2. —We argue that for a single disorder re¬ 
alization of (1) with 5 > 1, the eigenstates at larger 
(in absolute value) energy densities are more likely to be 
localized. We use the excited-state real-space renormal¬ 
ization group (RSRG-X) technique, which is a heuristic 
approach to solve or make progress towards solving the 
eigenvalues and eigenstates in strongly disordered sys¬ 
tems [12, 23, 30]. As any other renormalization group 
procedure, RSRG-X sequentially “integrates out” some 
degrees of freedom, while generating effective interactions 
between the remaining degrees of freedom. Thus, the 
disorder strength is modified. As the implementation of 
RSRG-X depends on the energy of the targeting eigen¬ 
state, the effective Hamiltonians for eigenstates at differ¬ 
ent energy densities may have different disorder strength. 

We do one layer of RSRG-X for (1). Assuming 
h » 1, we identify a set T C Z such that (i) \hi\ 
max{|/ij_i|, 1} for Vi € T; (ii) |T| is of the order 

of the system size. Fix 5=1. For each i £ T, (i) im¬ 
plies that the spin i is (approximately) in an eigenstate 
|sf = 0, ±1) of 5,f. Then a perturbative calculation shows 
that in the effective Hamiltonian acting on the remaining 
spins Z \ T, the terms in a neighborhood of i read 

- (/ii_i ^ s?)5?_i - (h i+1 - *l)5f., + H t _ M+1 , (3) 

where ||77i_i ji+ i|| = 0(l/\hi\) acts on the spins i — 1 and 
i + 1 so that the spin chain is not completely decoupled. 
Unlike sf = 0, sf = ±1 induces fluctuations in the ran¬ 
dom fields and thus increases the disorder strength in the 
effective Hamiltonian. We expect that I s ? I tends to 
be larger when the targeting eigenstate of (1) has larger 
(in absolute value) energy density; see the paragraph con¬ 
taining (13) for a quantitative analysis of a very similar 
statement. 

Since the disorder strength of the parent or effec¬ 
tive Hamiltonians for individual eigenstates depends on 
energy density, we propose “energy-resolved disorder 
strength” as a physical mechanism for the appearance of 
mobility edges. Obviously, our arguments apply to not 
only (1) but also a large class of models, which suggests 
the universality of this mechanism. We will construct a 
simple spin chain model (5) of MBL such that the energy- 
resolved disorder strength can be solved analytically: It 
decreases continuously and monotonically as energy den¬ 
sity increases. This implies a mobility edge below (above) 
which almost all eigenstates are localized (delocalized). 

Preliminaries.- To proceed, we need some formal def¬ 
initions. An (interacting) model is MBL if almost all its 
eigenstates are localized. In the literature, there are sev¬ 
eral (inequivalent) criteria for whether an individual state 
is localized. Our construction does not rely on which cri¬ 


terion to use. For example, you might keep in mind that 
a state is localized if it satisfies an area law for entangle¬ 
ment [5, 9, 26]. 

We take a detour and discuss QCLT on the distribution 
of the eigenvalues of a local Hamiltonian. In our context, 
it says that all but an exponentially small (in the system 
size) fraction of eigenstates have the same energy den¬ 
sity up to corrections vanishing in the thermodynamic 
limit. QCLT can be proved for any local Hamiltonian 
[6, 8, 15]. Here we just prove (a weak version of) it in all 
weakly interacting systems. Indeed, the single-particle 
spectrum of any (homogeneous or disordered) local free- 
fermion Hamiltonian is bounded. For a random rnany- 
body eigenstate, its energy is a random sum of single¬ 
particle energies, and the (classical) CLT implies that 
its energy density is close to the mean with overwhelm¬ 
ing probability, cf. ID random walk with random but 
bounded step size. (Weak) interactions can broaden any 
energy density interval by at most e if the norm of the 
interaction terms per site is upper bounded by e. 

The definition of MBL reads “... if almost all ...” be¬ 
cause one cannot rule out the possibility that a very small 
fraction of eigenstates are delocalized (especially in ran¬ 
dom systems). However, this seemingly innocent defini¬ 
tion has a caveat. Even if all but an exponentially small 
fraction of eigenstates are localized, there can still be an 
exponential number of delocalized eigenstates. Further¬ 
more, it is possible that almost all eigenstates away from 
the mean energy density are delocalized or that mobility 
edges exist. Indeed, QCLT implies that a model is MBL 
if almost all its eigenstates at the mean energy density 
are localized. 

As we will encounter probability distributions of var¬ 
ious shapes, a formal definition of disorder strength ap¬ 
pears necessary. Let A be a real-valued random variable 
with probability density function p(x), where p(x) > 0 
and J ^p(x)dx = 1. We only consider symmetric 
probability distributions with compact support such that 
p(x) = p(—x) and j^ A p(x)dx = 1 for A = 0(1). The dis¬ 
order strength s(X) of a random variable is a measure of 
how broad p{x) distributes. Formally, s is a nonnegative 
function of random variables satisfying 

(i) s(X) = 0 if and only if p(x) is a delta function. 

(ii) Let Xk(k = 1, 2) be a uniform random variable on the 
interval [— A*,, A*,] with Ai < A 2 , and (1 — A)XlOrAA 2 
denote a random variable which is X 2 with probability 
0 < A < 1 and X\ otherwise. Then s((l — A)XlOrAX 2 ) 
is a continuous monotonically increasing function of A. 

(iii) Other postulates irrelevant to us. 

A canonical example of s is the variance of a random 
variable, but this may not be a faithful measure of dis¬ 
order strength from a physical point of view, e.g., in the 
random-field Heisenberg chain (1), two different distri¬ 
butions of hi s with the same variance may correspond 
to the different phases. Hence, we do not specify s, but 
only assume that a faithful measure of disorder strength 
exists for the random spin model (4). 

The model. —We now present the details of our model. 
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Similar but not identical constructions appeared previ¬ 
ously in different contexts, e.g., [11, 22, 24, 28, 34]. 

We start with a disordered Hamiltonian on a spin chain 
(in the thermodynamic limit), where the local dimension 
d = 0(1) of each spin is a small constant, 

H = E Hij+i + XiHi. (4) 

Z 

Here, ||H^ i+1 || < 1 is a translationally invariant nearest- 
neighbor interaction between the spins i and i+ 1; ||.ffj|| < 
1 is a translationally invariant on-site term acting on the 
spin i; Ai’s are i.i.d. random variables. Without loss of 
generality, we assume tr = tr Hi = 0 so that the 

mean energy density of H is 0. 

Suppose H has a many-body localization- 
delocalization transition tuned by disorder strength: 
There is a critical s* = 0(1) such that almost all eigen¬ 
states of H are localized (delocalized) if s(Xi) > (<)s*. 
Let Xk(k = 1,2) be a uniform random variable on 
the interval [—A*,, A*.] with s(Ai) < s* < s(A 2 ). As 
s((l — A)AiORAA 2 ) is a continuous monotonically 
increasing function of A, there exists 0 < A* < 1 such 
that s* = s((l — A*)XiOrA*X 2 ). Although s* is fixed, 
A* is tunable in the sense that A 2 is tunable. 

Based on (4), we construct a simple spin chain model 
(5) of MBL with a tunable mobility edge. Here, a mobil¬ 
ity edge is an energy density below (above) which almost 
all eigenstates are localized (delocalized). We have two 
spins per unit cell (or two particle species in the language 
of [11, 22, 24, 34]). The main idea of the construction is to 
have the second spin species control (via the second term 
in (5)) the disorder strength in the effective Hamiltonian 
(7), which acts on the first spin species and depends on 
the state of the second species. Then, we move the eigen¬ 
states of (5) whose effective Hamiltonians (7) have strong 
(weak) disorder to the bottom (top) of the energy spec¬ 
trum of (5) by adding a uniform field (the third term in 
(5)) on the second spin species, which results in a nega¬ 
tive gradient in the energy-resolved disorder strength. 

In each unit cell, the first spin has local dimension d, 
and the second has local dimension 2. You may combine 
these two spins into a single spin of larger local dimension 
2d = 0(1) if you prefer one spin per unit cell. We label 
all spins by two indices (i. j). Here, i S Z is the unit cell 
index; j = 1,2 is the species index. The Hamiltonian is 

H = ^ + H j,i,j, 2 + of, (5) 

iez 


i - r 

•> 


r t 

'*•' 'tv 



FIG. 1. (Color online) We have two spins (dots) per unit 
cell (dashed circle). Conditioned that the right spins (second 
species) are in the state | • • • ftt-J-tt-ft • • •) (black arrows), 
the lengths and directions of the blue arrows, respectively, 
illustrate the magnitudes and signs of the random fields acting 
on the left spins (first species) in a typical disorder realization 
of the effective Hamiltonian (7). The short dashed lines show 
the maximum random-field strength in every unit cell, which 
is “controlled” by the state of the second spin species. 


As {of} is a set of conserved quantities, in any eigen¬ 
state of H the spins (* € Z, 2) are in a product state. For 
a particular configuration {erf } € {1,— 1}® Z , the effec¬ 
tive Hamiltonian acting on the first spin species reads 


H(K}) = E 


eH, 


i,i +1 


£ A( 3 - rT z)/ 2 H i 


Z 


1 + a 2 




( 7 ) 


which is equivalent to (4) after rescaling. Clearly, the 
second spin species controls the disorder strength in the 
sense that X( 3 _ a *y 2 = X ± if of = 1 and X( 3 _ a *y 2 = A 2 
otherwise. As of = ± 1 with equal probability in a ran¬ 
dom eigenstate of H, the (overall) disorder strength in 
(7) is s(^AiOR^A' 2 ). Compared with the critical disor¬ 
der strength s* = s((l — A*)A 1 ORA*A 2 ), H is MBL if we 
assume A* < 1/2. 

Let z := of/n be the average magnetization of the 
second spin species, where n is the number of unit cells. 
The energy density (per unit cell) of an eigenstate of H 
is well approximated by 2 . Indeed, 

|H({of»/n -*\<-£ l|H -'‘f( + Ail|H ‘ 11 < e . (8) 

n z ' 1 + A 2 

i 

Furthermore, QCLT says that for any particular config¬ 
uration {erf} almost all eigenstates of H({of}) have en¬ 
ergy density z. As of = ±1 with equal probability in 
a random eigenstate of H, 2 follows a rescaled Bernoulli 
distribution, and 1 2 1 = o(l) for almost all eigenstates of 
H. This allows the presence of a mobility edge at energy 
density e > 0. 

We now solve the energy-resolved disorder strength in 
our model by (approximately) identifying 2 with energy 
density. When 2 is fixed, 


where H ijyj* acts on the spins (i,j) and ( i',j'), and 
of is the spin-1/2 Pauli operator for (i, 2). Specifically, 
H *,i,*+i,i — e 7Lj j_)_i/( 1 + A 2 ) is translationally invariant, 
and Hj,iy, 2 = eHi ® diag{Ai, A 2 }/(1 + A 2 ) involves ran¬ 
domness, where e is a very small constant. (5) is a weakly 
interacting system as 


IIHjj^+ipll + ||Hj )1)ii2 || < 1 ^ — e ' 


X(3-trf)/2 — 2 AlOR 2 A 2 . (9) 

Hence s(A( 3 _ cr *)/ 2 ) is a continuous monotonically de¬ 
creasing function of 2 . Note that by fixing 2 the ran¬ 
dom variables A( 3 _ cr? w 2 ’s in different unit cells become 
slightly correlated, but we expect no significant physical 
effects of this very weak correlation. Compared with the 
critical disorder strength s* = s((l — A*)AiOrA*A 2 ), a 
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mobility edge appears at the energy density e such that 

^(-^(3—cr?)/2)| z=e ~ S* = * r 6 = 1 2A* > 0, (10) 

and almost all eigenstates of H below (above) the energy 
density e are localized (delocalized). 

Discussions & extensions. —We have established a sim¬ 
ple reduction from a random spin model (4) to another 
(5) with a nonzero gradient in the energy-resolved disor¬ 
der strength. If the former has a many-body localization- 
delocalization transition tuned by disorder strength, then 
the latter has a mobility edge. Furthermore, we expect 
that the transition tuned by disorder strength in the for¬ 
mer and that tuned by energy density in the latter have 
the same critical exponents. It should be clear that such a 
reduction can be established in general. Although we as¬ 
sumed that (4) is a (random) spin chain with only on-site 
disorder, a minor modification of our construction works 
in any spatial dimension and applies to cases where the 
original model has bond disorder. 

The model (5) has a mobility edge below (above) which 
almost all eigenstates are localized (delocalized). Adding 
a minus sign (trivially) reverses the energy spectrum so 
that almost all eigenstates above (below) the mobility 
edge are localized (delocalized). To construct a spin chain 
model of MBL with two mobility edges, it suffices to mod¬ 
ify only the third term in (5): 

H 2 = Y, H bM+1,1 + + hi(I + at )/2, (11) 

Z 


where hi = ±1 with equal probability. Note that {of} 
remains a set of conserved quantities, and the effective 
Hamiltonian acting on the first spin species becomes 

tr /r zn e -^>i+l + e ^(3-of)/2#i hi(l + of) 

l) = X-rnw- —2—■ 

ie z J 

( 12 ) 

Similarly, H 2 is MBL if A* < 1/2. 

Indeed, the effect of the third term in (11) is to move 
the eigenstates of (11) whose effective Hamiltonians (12) 
have weak disorder to the ends of the energy spectrum 
of (11), which results in a “bump” in the energy-resolved 
disorder strength. As before, we (approximately) identify 
y := /&»(! + erf)/(2n) with energy density. When y is 

fixed, the probability distribution function of z is 


P{z ) 


C 


2 Z «/ 2 


^ (l-z)re ^ | ^ (l+z+2y)n 


where C is a normalization factor, 
peaked at z = y 2 , (9) becomes 


^ i ^ ( 1 + z ~ 2 y) n j i ’ 

(13) 

As p(z) is sharply 


most all eigenstates of H 2 at energy densities (in abso¬ 
lute value) < (>)-\A — 2A* are localized (delocalized). It 
is also straightforward to construct a random spin chain 
model with two mobility edges such that one region of de¬ 
localized eigenstates sandwiches two regions of localized 
eigenstates. 

All our constructions can be made translationally in¬ 
variant. We show how to do this using (5) as an ex¬ 
ample. Recall that Af(/c = 1,2) is a continuous uni¬ 
form random variable on the interval [—Afc,Afe], Let 
X k be a discrete uniform random variable on the set 
Sk = {xk.i,Xk, 2 , ■ ■ ■ ,Xk,d'} such that (i) d! = 0(1) is 
a constant; (ii) \xkj\ < A k for j = 1,2,... ,d'; (iii) 
Sk = —Sk so that the distribution of X' k is symmetric; 
(iv) The numbers in Sk are generic in the sense that (16) 
has no accidental degeneracy; (v) s(Xk) = s(X’ k ) so that 

H' = ^H iiM+lil +H' ;Mi2 + uf (15) 
ie z 

reproduces the physics of (5), where H' t 2 = eH, (g) 
diag{A{, A'}/(1 + A 2 ). 

We encode all disorder realizations of (15) into a single 
translationally invariant Hamiltonian (16) by adding a 
third spin to every unit cell. These spins are of local 
dimension d ', and are labeled by (i, 3) for i £ Z. It 
suffices to modify only the second term in (15): 

H T = ^H. i , M+1>1 + H'+df, (16) 

Z 

where H' is a translationally invariant term acting on all 
three spins in the unit cell i. In the computational basis, 

H' = tHi® diag{xi a ,a:i j2 ,... ,a;i id -,a: 2) i,a: 2)2 ,.. .,x 2 ,d'}- 

(17) 

Let fj = diagjl, 2,..., d'} be an operator acting on 
the spin (i, 3). As {t/} is a set of conserved quan¬ 
tities, in any eigenstate of Hr the spins (i € Z, 3) 
are in a product state. For a particular configuration 
{t,} e {1,2,..., d'}® z , the effective Hamiltonian acting 
on the first and second spin species is a disorder realiza¬ 
tion of (15), where the random variable X' k in the unit 
cell i takes the value Xk, Ti - Thus, Hr is a translationally 
invariant Hamiltonian with a mobility edge. 

Outlook. An important open problem is whether the 
models we constructed are stable against perturbations. 
We now lack the tools to make progress on this problem. 
A proof of stability (even merely for a particular model) 
is notoriously difficult. Numerical methods (especially 
ED) are usually limited to small system sizes, and may 
suffer from strong finite-size effects [22]. 


ACKNOWLEDGMENTS 


1 + ir 1 — IT 

^■(3—<rf)/2 = --AfiOR--- X 2 . (14) 

with overwhelming probability. Compared with the crit¬ 
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